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Abstract. This paper deals with Student t-processes as studied in (Cufaro Petroni 
2007). We prove and extend some conjectures expressed by Cufaro Petroni about the 
asymptotical behavior of a Student t-process and the expansion of its density. First, 
the explicit asymptotic behavior of any real positive convolution power of a Student 
t-density with any real positive degrees of freedom is given in the multivariate case; 
then the integer convolution power of a Student t-distribution with odd degrees of 
freedom is shown to be a convex combination of Student t-densities with odd degrees 
of freedom. At last, we show that this result does not extend to the case of non-integer 
convolution powers. 
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1. Introduction 



In a recent contribution (Cufaro Petroni 2007), N. Cufaro Petroni derived several results 
about the behavior of some non stable Levy processes with Student t-marginals. More 
precisely, he considered the random walk 

N 

Z N = X,i 

i=l 

where N G N and each independent step Xi follows a Student t-distribution with 
/ = In + 1 degrees of freedom, n G N. We recall that the Student t-density with 
f = 2v> degrees of freedom [y > 0) is 

/„(*) = Mi + *T<"»; ^ = ^0y 

The family of Student t-densities includes the Cauchy density for / = 1 and the 
scaled density -^=f u (-^=) converges to the Gaussian density as / — > +oo. All Student 
t-distributions are heavy tailed. Grosswald (Grosswald 1976) proved that they are 
infinitely divisible. They also have the stronger property of being self-decomposable, cf. 
(Steutel & van Harn 2004). 

Stochastic processes with Student t-marginals and various types of dependence 
structures have been proposed in (Heyde & Leonenko 2005), most of them with 
dependent increments. On the other side, Cufaro Petroni's paper deals with Levy 
Student t-processes, which exist by the infinite divisiblity of the Student t-distribution. 
In both cases, these processes have applications in finance (Schoutens 2003) and in 
physics (Vivoli et al. 2006). 



2. Three conjectures by Cufaro Petroni 

Cufaro Petroni obtained precise results about the process only in the case of / = 3 
degrees of freedom; however, he expressed three conjectures about the extension of these 
results to more general cases; the first conjecture is 

Conjecture 1. For all N G N and for all f = 2n + l, n G N, the distribution of N~ l Z^ 
is a convex combination of Student t-distributions with odd degrees of freedom. 

The two remaining conjectures concern the distribution of Z N for non-integer values 
of N, which makes sense because the Student t-distribution is infinitely divisible: the 
c— fold convolution of the distribution f v is defined, for any real positive c, as the inverse 
Fourier transform 

1 / >+ °° • 
f?( x ) = to J e-[Mu)] c du, 

where ip u (u) is the characteristic function of the Student t-distribution 

ip v (u) = k u (H) 
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with 

K (u) = Y^f K ^ (m) ' u > °- (1) 

Here K v is the modified Bessel function of the second kind also called the Macdonald 
function. The expression ([T]) reduces to elementary functions exactly when v = 
n + 1/2, n — 0, 1, . . . because 

k n+ i(u) = e- u q n (u),u > 0, (2) 

where q n is a polynomial of degree n with positive coefficients, called the n'th Bessel 
polynomial. It is given as 



n 



q n (u) = yy^u k , (3) 

where 



k=0 



(n) _ (-n) k 2 k 

k (-2n) k k\- {i) 



The first examples of these polynomials are 



u 2 



q (u) = 1, qx (u) = 1 + u, q 2 («) = 1 + U + y, 

cf. (Berg & Vignat 2008). 

Cufaro Petroni's second conjecture concerns the asymptotic behavior of the density 
of the c— fold convolution of f v : 

Conjecture 2. For every c > and v > 0, i/ie asymptotic behavior of the c—fold 
convolution f* c is given by 

cA 

Cufaro Petroni's last conjecture is an extension of Conjecture 1 to the c—fold 
convolution f* c as follows: 

Conjecture 3. Conjecture^ extends to non-integer sampling times c under the following 
form: for all u > and all c> 0, 

f +OQ 1 x 

J i/ C C 

/or some distribution q VQjC (i/). 

In this paper, we show that conjecture ([1]) holds true and give an extended version of 
it; likewise, we prove the conjecture (j2J). We were unable to prove or disprove conjecture 
([3]), but we disprove a discrete version of it in the case where v = n + \ with n G N. 

Moreover, we consider in the rest of this paper the multivariate context: all Student 
t-variables are supposed rotation invariant d— dimensional vectors. The multivariate 
Student t-density is given, for x = (xi, . . . , x d ) e K d by 

/„(x) = A d , p (1 + | x | 2 )-^ /2 , A d , - J^±SL 
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d 

(x,y) = ^2xiy t , |x| = (x,x)5, x,yGR d . 
i=i 

3. First conjecture: the odd degrees of freedom case 

Cufaro Petroni's first conjecture (Cufaro Petroni 2007, Prop. 5.2) is that if Xi is a set 
of independent Student t-distributed random variables with / = In + l,n G N degrees 
of freedom, then the density of the distribution of the normalized iV— th step of the 
random walk 

1 - 

i=i 

writes as 

nN 
k=n 

with (5^' N ^ > 0, n < k < nN. We extend and prove this conjecture as follows: 
Theorem 1. If N G N and 

N 

i=i 

where are positive numbers with sum 1 and Xj are independent d—variate Student 
t-distributed, each with f\ = 2rij + 1 (n^ G N) degrees of freedom, then the density ofY N 
is 

ni+...+n N 
j=mm(ni,...,njv) 

where the coefficients /3j are nonnegative with sum 1 and do not depend on the dimension 
d. 

Proof. The characteristic function of the d— variate Student t-distribution is 
<p(u) = fc„(M) 

where the function k u is given by ([T|). Since for z/j = rij + |, this function reads 
fc^dul) = e _ l u l9 ni (|u|) where g ni is the Bessel polynomial of degree n u the result follows 
from (Berg & Vignat 2008, Th. 2.6): 

L 

q ni (aiu)q n2 (a 2 u) ■ ■■q nN { a Nu) = ^fyq^u), u G R (6) 

j=i 

with nonnegative coefficients (3j with sum 1 and I = min(ni, . . . , tin), L = n% + ■ ■ ■ + 
n N . □ 
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As a particular case, choosing aj = i, 1 < i < N for AT G N, we deduce that the 
coefficients ^ in (j5J) are positive, and thus the density <? is a convex combination of 
Student t-distributions with odd degrees of freedom. 

We are not able to provide an expression for the coefficients (5^ ,N ^ which can be used 
directly to see the non- negativity. Using Carlitz' formula, see (Berg & Vignat 2008), 

n 

u n = Y t tf? ) q j {u) t n = 0,l,... (7) 

3=0 



with 



| (n+l)! (-l)-'(2j)l f n-i < • < 

= J 2" (n-i)!i!(2i+l-n)! 101 2 — — (g) 

1 for < J < afl 



it is possible to write 



N m n N / N \ fciH hfcjv 

iim«v)=e-e nw e *." =,+ +tal 

i=i fci=o fc]v=o \i=i / i=o 



Qi[u) 

(n) 



which gives an expression for (3j in (joj), but because of the varying sign of 5j , it is not 
possible to see directly that f3j > 0. 

If a,j = 1/N and n\ = . . . = tin — 1, i-e. the case of / = 3 degrees of freedom where 
qiiu) — 1 + u, this formula simplifies to the expression given in (Cufaro Petroni 2007, 
Prop. 5.2). It is claimed that the expression is positive, but no convincing argument is 
given. 

4. Second conjecture: the asymptotic behavior of the Student process 

A second property studied by Cufaro Petroni is the asymptotic behavior of the 
distribution of the random walk Zj^; in the case of / = 3 degrees of freedom [y = |) , 
he obtains the following result (Cufaro Petroni 2007, Prop. 5.1): for all c > 0, 

2c 

flfr) 7, X -> +00. 

We provide now an extension of this result to any value / = 2u (i> > 0) of degrees of 
freedom. Cufaro Petroni's argument is via Fourier analysis. This argument becomes 
very technical if one tries to generalize it to arbitrary degrees of freedom. Our proof is 
based on results about subexponential distributions. 

Theorem 2. For any c > and any v > 0, the density of the c—fold convolution of the 
d—variate Student t- distribution behaves asymptotically as 

rA , 

£*C I \ ^^d,!/ I I 

f» ( x ) ~ l x l -> +°°- 

Proof. The proof is based on a series of lemmas given in the last section. The d— variate 
Student t-distribution is subordinated to the d— variate Gaussian semigroup 



<fe(x) = (4vrt)- d / 2 exp (~^\ t > 0, x G 
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by the inverse Gamma density, i.e. 

/„(x)= / g t (x)dH v (t), 
Jo 

where H v (t) is the inverse Gamma distribution with density 
K (*) = C u exp {-\- ) r"- 1 , t > 0, a ' 



4t y " 2 2 -r(z/)' 

From this representation we deduce in Lemma [1] the same representation for the c-fold 
convolution power of the Student t-density, namely 

" + 0O 



f? (x) = / g t (x) dH? (t) . 

JO 



We note that this property is very general, in the sense that it holds for any infinitely 
divisible probability distribution dH (t) on [0, oof . 

The next step of the proof is the derivation of the asymptotic behavior of the c-fold 
convolution power of the inverse Gamma density h u (t) : by Lemma [2j this reads 

hl c (t) ~cC v t~ v - x , t->+oo. 

Finally, we show in Lemma H] that this asymptotic behavior implies, by 
subordination to the Gaussian semigroup, the desired asymptotic behavior of the c- 
fold Student t-convolution. □ 

As a consequence of this theorem, we deduce the following 

Corollary 1. In the case where the number of degrees of freedom 2v = 2n + 1 is an odd 
integer and with integer N, the coefficient (3^' N ^ in (FJP reads 

1 



(n,JV) 



AT 2 ™ 



Proof. Since the coefficients (3^^ do not depend on the dimension d, we consider the 
case d = 1. The function v in (j3J) is the density of the normalized random walk i J2f=i Xi 
and thus writes 



v (x) = Nf*"i (Nx) . 



l 

2 

By Theorem [21 



v (x) ~ N 2 K+ l +2 = ^±ix- 2 ™- 2 , x - +oo. 

Since each Student t-distribution f k+ i in (jSJ) has asymptotic behavior 

f k+ i (x) ~ A k+ ,x- 2k -\ 
we deduce that 

v (x) ~ (3^A nH x~ 2n - 2 . 
Identification of the two equivalents yields the result. □ 
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5. Third conjecture: non integer sampling time and odd degrees of freedom 

In this section, we prove by contradiction the following result 

Theorem 3. For all c > 0, c ^ N and v = n+ ~, n G N, the 1—variate density f* c can 
not be expanded as 

+oo 1 

i=o 

TOi/j parameters (3j > 0. 

Proof. We remark that integrating equality ([9]) over R yields X]j=o = 1 so th & t the 
sequence {(3k) is summable. The Fourier transform of (jUJ) reads 

+ 0O 

fe n+±H = ^ fa exp(-cu)qj(cu), u>0 

j=0 

where qj is the Bessel polynomial of degree j. Thus, by LemmaEl the sum Y^j^o PjQ.j( cu ) 
is an entire function, so that the function k c , 1 (u) extends to an entire function. But 

K+d u ) = exp(-cu) [q n (u)} c , 

and since c is not an integer, the function q^ is not holomorphic at any of the complex 
roots of q n , what concludes the proof. □ 

6. Lemmas for the proof of theorems [2] and [3] 

Lemma 1. The c—fold convolution of the density f v reads, for all c> and v > 0, 

r+oo 

ft c (x) = / g t (x) dH? (t) . 
Jo 

Proof. Let us consider the d— dimensional Fourier transform T of the right hand side 



gt (x) dH? (t) 



o 



(y) 



rOO 



F[g t {*)]dHl c {t) 



o 



(y) 



r+oo 

/ exp (-t\y\ 2 ) dH*/ (t) . 
Jo 



This last integral is nothing but the Laplace transform £ of H* c evaluated at |y| 2 , and 
thus coincides with 

(r+oo \ c 

J exp(-t| y | 2 )^(t)J =(^(/,(x))) c ( y ). 

The result follows by considering the inverse Fourier transform. □ 

Lemma 2. For all c > and v > 0, i/ie c—fold convolution of the inverse Gamma 
density has asymptotical behavior 

hl c (t) ~cC u r v ~\ t^+oo. 
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k (<) = aexp f-^j t-"- 1 ~ at"" -1 , * -> +00, 

the tail function H v (t) = 1 — H v (t) of the inverse Gamma distribution has the asymptotic 
behavior 

S v it) ~ -r", t -> +00. 

This tail function is thus regularly varying and, by (Feller 1950, p. 278), 
Hf ( t ) „ ?Q. t -» 



where H* 2 is the tail of H* 2 , so that the inverse Gamma distribution is subexponential in 
the sense of Chistyakov, cf. (Chistyakov 1961). Since it is moreover infinitely divisible, 
we deduce by (Embrechts et al. 1979, corollary 1, p. 340) that for all c > 

rC 

Hi c (t)~—r v , t^+00. 

By Lemma [3] below the density h* c (t) is ultimately decreasing, so the result follows by 
application of the monotone density theorem (Bingham et al. 1987). □ 

Lemma 3. The c-fold convolution of the inverse Gamma density is ultimately 
decreasing. 

Proof. The inverse Gamma density h v (t) is a generalized Gamma convolution, and so is 
the convolution power /i* c , cf. (Steutel & van Harn 2004, p. 350). Since its left extremity 
is 0, we deduce from (Steutel & van Harn 2004, prop. 5.5) that it is unimodal, and thus 
ultimately decreasing. □ 

Lemma 4. The asymptotic behavior of the c-fold Student t-convolution is 

- W ~ j^p5> l x l -> °°- 
Proof. Since by Lemma [2] 

K c (t) ~ cC u r v -\ 

for any a and b such that a < cC u < b, there exists to > such that for all t > to? 
at-"- 1 < h* v c {t) < bt-"- 1 . 

From 

+OO pto r+00 

g t (x) C (t) dt= g t (x) h* v c (t) dt + g t (x) h? (t) dt, 

JO Jt 
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it follows that 



to f + OO f+OO 



Ul fx) /,;; (/) <it 4 / gt ( x )—dt<J^ g t (x)hl c (t)dt 



< ^ ft (x) h*; (t) dt + jf & (x) —dt. 
But the integral C° g t (x) /i* c (t) dt is o(|x| _2l/_d ) for |x| — > oo because 



-to 
'0 

jT ft (x) c (o dt = jT ^-Jp ex P (t) dt 

|2\ /-to ]_ 



A simple computation gives 

"+00 r)2i^ 



/ g t hc)t- l/ - 1 dt= \x\- 2v - d ^ / exp (-wK+^dw 



't « ' JO 



hence 



/"+00 /-+0O L 

limsup |x| 2l/+d / ft (x) (t) dt < limsup |x| 2l/+d / ft (x) — dt 

|x|^oo JO |x|^oo Jt " 



b2 2l/ ( d 
-T ( zy + - 



d/2 V 1 2 



and 



7T 



/•+oo p+oo r 

liminf |x| 2l/+d / ft (x)/i* c (t)dt > liminf |x| 2l/+d / ft' ' 

|x|-*oo J ' |x|-»oo y tQ 



x) — — dt 

I t u+l 



«2 2t ' / d 
-r ( zy + - 



7T 



d/2 ' 



so that finally 



J ft (x) k c (t) dt ~ c -^r lu + - a|x| 



Hi/ 1 

□ 

Lemma 5. For fixed k, the coefficients a k of the Bessel polynomial of degree n > k is 
increasing in n and 



lim af ] = — . 

n-*+oo K kl 

Proof. From (j3j) we get 



1 



fc-i 

(n) 1 TT n ~ J ^ 1 



'■• ~ fc!ll„_2 - fc!' 

j=l 2 
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where each of the (k — 1) terms of the product 



n ~ J _ i 2 



is increasing and converges to 1 with n. □ 
Lemma 6. Consider the infinite series (S) equal to 

+0O 

k=0 

where qk{z) is the k 'th Bessel polynomial and c^gC. Then the three following assertions 
are equivalent: 

(i) (S) is absolutely convergent for z = 0, 

(ii) the sequence (c^) is absolutely summable, 

(Hi) the series (S) converges absolutely and uniformly on any compact subset of C. 

Proof, (i) =>- (ii) since q k (0) = 1. 
(ii) =>- (iii) since 



n iii. +oo I I u 

\c n q n (z)\ < \c n \q n (\z\) < \c n \^2^- < l c nl$Z^T - l c «l ex P(^) 

fc=0 ' k=0 

for some constant K. The first inequality holds since all the coefficients aj^ of the Bessel 
polynomial q n are positive; the second inequality is a consequence of the majorization 
— h P rove d in Lemma EJ the third inequality is straightforward, and the last 
inequality ensues from the assumption that z belongs to a compact subset of C. Since the 
sequence (c„) is assumed absolutely summable, the absolute and uniform convergence 
of (S) is a direct consequence of the above majorization. 

(iii) =^ (i) trivially. □ 
7. Conclusion 

In this paper, we have proved two of the conjectures as expressed by Cufaro Petroni in 
(Cufaro Petroni 2007), and disproved a simple version of a third conjecture. We note 
that these results extend naturally to d— dimensional Student t-vectors with correlated 
components, with density 

/,(x) = ^(l+x^- 1 x)-^f) 

and characteristic function 

9?(u) = k v (V u*i^u) 

where K is a symmetric and positive definite matrix; Theorem [1] holds unchanged and 
the asymptotic result of Theorem [2] still holds by replacing constant A djU by — 
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